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Abstract 


High transverse momentum jets produced in pp collisions at a centre of mass energy of 7 TeV 
are used to measure the transverse energy-energy correlation function and its associated 
azimuthal asymmetry. The data were recorded with the ATLAS detector at the LHC in the 
year 2011 and correspond to an integrated luminosity of 158 pb”L The selection criteria 
demand the average transverse momentum of the two leading jets in an event to be larger 
than 250 GeV. The data at detector level are well described by Monte Carlo event generators. 
They are unfolded to the particle level and compared with theoretical calculations at next- 
to-leading-order accuracy. The agreement between data and theory is good and provides a 
precision test of perturbative Quantum Chromodynamics at large momentum transfers. From 
this comparison, the strong coupling constant given at the Z boson mass is determined to be 


cusimz) - 0.1173 + 0.0010 (exp.) ;!!ooo 26 (tbeo.). 
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1 Introduction 


The study of jet production at the LHC provides a quantitative test of Quantum Chromodynamics, QCD, 
at the highest momentum transfers. Theoretical calculations for jet cross-sections in hadronic collisions 
have been carried out up to next-to-leading order (NLO) accuracy in the strong coupling constant as [1-3] 
and extensively compared with the data [4-10]. These calculations are valid for configurations with up to 
four jets in the final sfafe. 

Evenf shape variables have been measured in all major e'^e~ experimenfs, as well as in experimenfs 
af fhe elecfron-profon collider HERA. These sfudies were recenfly exfended fo hadron colliders wifh 
measuremenfs of fhe fransverse fhrusf and fhe fransverse minor [11, 12] af fhe Tevafron [13] and fhe EHC 
[14, 15]. 

Energy-energy correlafions (EEC), i.e. measuremenfs of fhe energy-weighfed angular disfribufions of 
hadron pairs produced in e'''e~ annihilafion, were proposed in Refs. [16, 17] as an alfernafive evenf shape 
variable nof based on fhe deferminafion of fhe fhrusf principal axis [18] or fhe sphericify fensor [19]. 
The EEC funcfion and ifs asymmefry, AEEC, were subsequenfly calculafed in 0(a^) [20-22], and fheir 
measuremenfs [23-35] have had significanf impacf on fhe precision fesfs of perfurbafive QCD and in fhe 
deferminafion of fhe sfrong coupling consfanf in e'^'e^ annihilafion experimenfs; a recenf review is given 
in Ref. [36]. The EEC are by consfrucfion nof affecfed by soff divergences, and as a consequence of fhis 
fhey are calculable af high orders. 

The frans verse energy-energy correlafion funcfion, TEEC, and ifs asymmefry, ATEEC, were proposed as 
fhe analogous variables af hadron collider experimenfs in Ref. [37], where predictions fo leading order 
(EO) were also presenfed. The NEO correcfions were calculafed recenfly in Ref. [38] using NLOJet-i-i- 
[2, 3]. These calculafions allow for a numerical deferminafion of fhe NEO predictions for fhe TEEC and 
ATEEC, i.e. fhe coefficienfs of fhe second order polynomials in fhe sfrong coupling consfanf. They are 
used in fhis paper for quanfifafive precision fesfs of QCD including a deferminafion of fhe sfrong coupling 
consfanf. The TEEC is defined as: 

1 dl _ 1 r do¬ 

er d(cos(;6) (T ^ J dVTidXT;d(COS 0) 


where fhe sum runs over all pairs of jefs in fhe final sfafe wifh azimufhal* angular difference (p = A^,y 
and XT, = Eli /Ej is the transverse energy carried by jet i in units of the sum of jet transverse energies 
Et = Et;. In order to cancel uncertainties that are constant over cos 0 € [-1,1], it is useful to define 
fhe azimufhal asymmefry of fhe TEEC (ATEEC) as 

1 dZ^'^y” ^ 1 dl _ 1 dl 
o-d(cos^) o-d(cos^) ^ crd(cos^) 

This Eeffer presenfs a measuremenf of fhe TEEC and ifs associafed asymmefry using high-energy jefs. 


( 2 ) 


2 The ATLAS detector 

The ATEAS defeefor [39] is a mulfi-purpose parficle physics defeefor wifh a forward-backward symmef- 
ric cylindrical geomefry and a solid angle coverage of almosf 47r. 

* ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the i/-axis points 
upward. Cylindrical coordinates (r, ip) are used in the transverse plane, ip being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 0 as ?/ = - In tan(0/2). 
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The inner tracking system covers the pseudorapidity range \t]\ < 2.5, and consists of a silicon pixel de¬ 
tector, a silicon microstrip detector, and, for \r]\ < 2.0, a transition radiation tracker. It is surrounded by a 
thin superconducting solenoid providing a 2 T magnetic field along the beam direction. A high-granularity 
liquid-argon sampling electromagnetic calorimeter covers the region I 77 I < 3.2. An iron/scintillator tile 
hadronic calorimeter provides coverage in the range \t]\ < 1.7. The endcap and forward regions, span¬ 
ning 1.5 < |? 7 l < 4.9, are instrumented with liquid-argon calorimeters for electromagnetic and hadronic 
measurements. The muon spectrometer surrounds the calorimeters. It consists of three large air-core 
superconducting toroid systems and separate trigger and high-precision tracking chambers providing ac¬ 
curate muon tracking for I 77 I < 2.7. 

The trigger system [40] has three consecutive levels: level 1 (LI), level 2 (L2) and the event filter (EF). 
The LI triggers are hardware-based and use coarse detector information to identify regions of interest, 
whereas the L2 triggers are software-based and perform a fast online data reconstruction. Finally, the EF 
uses reconstruction algorithms similar to the offline versions wifh fhe full defecfor granularify. 


3 Monte Carlo samples 

Mulfi-jef production in pp collisions is represenfed by fhe convolution of fhe production cross-secfions 
for parfon-parfon scaffering wifh fhe parfon disfribufion functions. Monfe Carlo (MC) generators differ 
in fhe approximafions used fo calculate fhe underlying shorf-disfance QCD process, in fhe way parfon 
showers are builf fo fake info accounf higher-order effecfs and in fhe fragmenfafion scheme responsible 
for long-disfance effecfs. For fhis analysis, fwo differenl MC approaches are used, depending on whefher 
fhe underlying hard process is considered fo be 2 ^ 2 or mulli-legged. The generated evenfs are fhen 
processed wifh fhe ATLAS full defector simulafion [41] based on Geant4 [42]. 

The baseline MC samples are generated using Pythia 6.423 [43] wifh fhe mafrix elemenfs for fhe under¬ 
lying 2^2 processes calculated af FO using fhe MRST2007LO* parfon disfribufion funcfions (PDF) 
[44] and mafched fo fransverse-momenfum-ordered parfon showers. The AUET2B fune [45, 46] is used 
fo model fhe underlying even! (UE) and fhe hadronisafion follows fhe Lund siring model [47]. 

Additional samples are generated wifh Herwig-i-i- 2.5.1 [48], using fhe CTEQ 6.6 PDF [49] and fhe 
UE7000 fune for fhe underlying even! [50]. HERwio-r-r uses angular-ordered parfon showers, a clusler 
hadronisafion scheme and ils own underlying-evenl paramelerisafion given by Jimmy [51]. 

A differenl approach fo simulate mulli-jel final slafes is followed by Alpgen [52]. This approach is based 
on FO malrix-elemenl calculations for 2 ^ n mulli-parlon final slafes, wifh n < 6, inferfaced wifh Her- 
wiG-i-JiMMY [53, 51] fo provide fhe parfon shower, hadronisafion and underlying-evenf models. Alpgen is 
known fo provide a good descripfion of fhe mulfi-jef final slafes as measured by ATLAS [54]. 


4 Event selection and jet calibration 

The dala used in fhis analysis were recorded in 2011 slI ^[s - 1 TeV and collecled using a single-jel 
Irigger. If requires al leasl one jel, reconslrucled wifh fhe anli-k^ algorilhm [55] wifh radius parameler 
R = 0.4 as implemenfed in FastJet [56]. The jel Iransverse energy, Ej - E sin 6 , is required fo be grealer 
lhan 135 GeV af fhe Irigger level. This frigger is fully efficienl al reconslrucled Iransverse energies above 
240 GeV. Taking info accounf fhe prescale factor of fhis Irigger, fhe dala collecled correspond to an ef¬ 
fective integrated luminosily of Xeff - 158 pb“' [57]. 
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Events are required to have at least one primary vertex, with five or more associated tracks with trans¬ 
verse momentum pj > 400 MeV. If there is more than one primary vertex, the vertex maximising ^ is 
chosen. MC simulated events are subject to a reweighting algorithm in order to match the average number 
of interactions per bunch-crossing observed in the data. 

In the analysis, jets are reconstructed with the same algorithm as used in the trigger, the anti-k, algorithm 
with radius parameter R = 0.4. The input objects to the jet algorithm are topological clusters of energy 
deposits in the calorimeters [58]. The baseline calibration for these clusters corrects their energy using 
local hadronic calibration [59, 60]. The four-momentum of an uncalibrated jet is defined as the sum of the 
four-momenta of its constituent clusters, which are considered massless. The resulting jets are massive. 
However, the effect of this mass is marginal for jets in the kinematic range considered in this paper. 

The jet calibration procedure includes energy corrections for multiple pp interactions in the same or 
neighbouring bunch crossings, termed “pileup” in the following, as well as angular corrections to ensure 
that the jet originates from the primary vertex. Effects due to energy losses in inactive material, shower 
leakage, the magnetic field, as well as inefficiencies in energy clustering and jet reconstruction, are taken 
into account. This is done using an MC-based correction, in bins of rj and pj, derived from the relation 
of the reconstructed jet energy to the energy of the corresponding hadron-level jet, not including muons 
or non-interacting particles. In a final step, an in situ calibration corrects for residual differences in the jet 
response between the MC simulation and the data using momentum-balance techniques for dijet, y-i-jet, 
Z-i-jet and multi-jet final states. This so-called jet energy scale (JES) [61] is subject to uncertainties in¬ 
cluding those affecting the energy of well-measured objects, like Z bosons and photons. The total JES 
uncertainty is given by a set of independent sources, correlated in pj. The uncertainty in the pj of indi¬ 
vidual jets due to the JES increases from for I 77 I < 1.8, to 5% for 1.8 < It/I < 4.5. 

The selected events must have at least two jets with transverse momentum pj > 50 GeV and pseudorapid¬ 
ity \r]\ < 2.5. The two leading jets are further required to fulfil pxi + Pti > 500 GeV. In addition, jets are 
required to satisfy quality criteria that reject beam-induced backgrounds [62] , as well as criteria for the 
fraction of the momentum of tracks within the jet which arise from the primary interaction vertex. The 
number of selected events in data is 3.8 x 10^, with an average jet multiplicity (Vjet) = 2.6. The resulting 
distribution for (pTi + Pt 2)/2 extends up to 1.3 TeV with an average value of 305 GeV. 


5 Results at the detector level 


The selected events are used to measure the TEEC and its associated asymmetry ATEEC, as defined in 
Equations (1) and (2). The TEEC distribution for a sample of N events is obtained by calculating the 
cosines of the angles in the transverse plane between all possible pairs of jets in each event. Every pair 
(/, j) represents an entry in the distribution, which is then weighted with the normalised product of the 
transverse energies. The weights Wij are defined as 


Wij = XT,-XT; ^ 


EjiE'ij 

iI.kETkf 


(3) 


such that for a given event their sum is always unity, as the self correlations i = j are also taken into 
account. The resulting distribution is then divided by the number of events, which normalises it to unit 
area. This weighting procedure reduces the sensitivity to the jet energy scale and resolution. 

Eigure 1 shows the TEEC and ATEEC distributions along with comparisons to detector-level Pythia, 
Herwig and Alpgen expectations. The TEEC exhibits peaks at cos cf) - \ (self correlations) and near 
COS0 = -1, with a rather flat central region around cos0 = 0. These features are similar to those 
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observed in e^e~ annihilation, as described in Ref. [31]. The central region is expected to be dominated 
by hard radiation processes while multiple soft radiation is expected to be important in the cos (/>-+! 
regions. 

The description of the TEEC is good in the back-to-back region cos 0 - -1 for both Pythia 6 and Alpgen. 
Differences up to 10% are observed in the central part, while the region of small angles shows differences 
as large as abouf 15%. The description by Herwig-i--i- is poorer. The ATEEC exhibifs a sleep fall-off, 
which is reproduced by bolh Pythia 6 and Alpgen. Hf,rwtg-p-p shows some discrepancies as large as 
30%. 



cos (|) cos (|) 


Figure 1; The detector-level distributions for the transverse energy-energy correlation TEEC (left) and its asym¬ 
metry ATEEC (right) along with comparisons to MC model expectations. The uncertainties shown are statistical 
only. The first bin of the ATEEC distribution has a negative value and is therefore not included in the figure. 


6 Correction to particle level 

The dala are correcled lo Ihe parlicle level in order lo lake info accounl defector efficiencies and resol¬ 
utions. This allows a direcl comparison wilh Iheorelical calculalions, as well as wilh measuremenls of 
ofher experimenls. 

Particle-level jels are reconslrucfed in MC evenfs using all particles wilh average lifetime r > 10“^^s, 
including muons and neulrinos. The kinematic seleclion crileria are fhe same as for the detector-level 
distribution. The unfolding relies on a bin-by-bin correction given by the ratios of the particle-level to 
detector-level distributions in the Pythia AUET2B sample, which is then applied to the detector-level dis¬ 
tributions in data. To check the effect of bin migrations on the unfolding procedure, an iterative Bayesian 
method [63] as implemented in RooUnfold [64] is also used. The convergence criteria is fulfilled when 
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the linear sum over all bins of the absolute relative differences from one iteration to the next drops below 
10“^. The method converges after five iferafions. The differences befween fhe fwo approaches are negli¬ 
gible, compared fo fhe sfafisfical uncerfainfies, in fhe full range of cos (p. This is expecfed due fo fhe high 
azimufhal resolution of fhe jef axis, which is 10 mrad. 

The following experimenfal sources of uncertainly are considered for Ihis measuremenl: 

• Jet Energy Scale: The uncertainly due fo fhe jef energy scale (JES) [61] is calculaled using MC 
techniques by varying each jef energy and momenlum by one slandard devialion for each of fhe 
63 independenl sources of fhe JES uncerfainly, and propagated fo fhe TEEC. These uncerfainfies 
depend on fhe jef Iransverse momenlum and pseudorapidily. The folal uncerfainly due fo fhe JES 
is calculaled as fhe sum in quadralure of all independenl uncerfainfies. In order fo investigate fhe 
effecl of possible correlations befween JES sources in fhe analysis, fwo alternative scenarios wifh 
weaker and sfronger correlalions have been considered [61]. The impacl of fhe change of correla¬ 
tion configurations, as well as of fhe number of JES independenl sources, on fhe value of a^{mz) 
and ils experimenfal error is found fo be negligible. 

The values of fhe JES uncerfainly are typically asymmefric for bofh fhe TEEC and ATEEC disfri- 
bulions, alfhough fhe values for fhis asymmelry are small. Thus, fhe positive and negative parls of 
the uncertainty are independently summed in quadrature. 

The TEEC distribution has a total uncertainty of up to 3.5% from the JES sources, the largest 
contributions being due to close-by jets and to the different response to jets initiated by gluons or 
quarks. This is the dominant experimental systematic uncertainty in the analysis. 

• Jet Energy Resolution: The uncertainty in the jet energy resolution [65] is propagated to the 
TEEC by smearing each jet transverse momentum by a pj- and ? 7 -dependent factor accounting for 
the resolution uncertainty. The size of this uncertainty is below 1% for both the TEEC and the 
ATEEC distributions. 

• Pileup: The pileup uncertainty is estimated by comparing the ratio of the detector-level TEEC 
and ATEEC distributions obtained in samples with reduced (ju < 6) and enhanced pileup activity 
(ji > 6). Here jj. is the average number of interactions per bunch crossing [57]. These ratios are 
formed in both data and MC simulation and the difference is assigned as fhe pileup sysfemalic 
uncertainly, which is as large as 2% (4%) for fhe TEEC (ATEEC). The size of fhis dedicated 
esfimale is larger lhan whal is predicled by fhe sum of fhe fwo sources of uncerfainly due fo pileup 
included in fhe JES uncerfainly. The envelope of fhe fwo differenl esfimales is used. 

• Parton shower modelling: To estimate fhe uncerfainly due fo fhe parlon shower modelling, fhe 
dafa unfolded wifh Pythia 6 and Herwig-i--i- are compared. The parlon shower and hadronisalion 
models in fhe fwo generators are differenl, as is fhe implemenlafion of UE effecls. The size of fhis 
uncerfainly is as large as 3.5% (2.5%) for fhe TEEC (ATEEC). 

• Unfolding: To estimate fhe uncertainly associaled wifh fhe unfolding procedure, a dala-driven 
mefhod is used fo lesf ifs sfabilily. This melhod relies on fhe reweighfing of fhe parlicle-level pro¬ 
jection of fhe unfolding Iransfer malrix so lhal fhe agreemenl befween fhe deleclor-level projection 
and fhe dafa is enhanced. This modified deleclor-level dislribufion is Ihen unfolded using fhe cor- 
recfion facfors described above. The difference befween fhe modified parlicle-level dislribufion and 
the nominal one is then taken as the uncertainty. This uncertainty is smaller than 0.5% for the full 
cos (p range. 
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Other possible sources of uncertainty are also studied, such as the jet angular resolution and jet quality 
selection procedure. They are found to be at the per mille level, much smaller than the statistical uncer¬ 
tainty on the corrected data, and are therefore neglected. To reduce the effect of statistical fluctuations, 
all the independent systematic uncertainties discussed here are smoothed separately. Figure 2 shows the 
breakdown of the systematic uncertainties for both the TEEC and the ATEEC, together with the total, 
obtained as the sum in quadrature of every independent source discussed above. 




cos (|) cos (j) 

Figure 2: Relative systematic uncertainties for the TEEC (left) and the ATEEC (right) as a function of cos (p. 

The TEEC and ATEEC distributions, once corrected for detector effects, are shown in Eigure 3, together 
with their total uncertainties, while numerical values are given in Tables 1 and 2. 

As already seen in the detector-level distributions, Pythia 6 and Alpgen give a fair description of the 

data both for the TEEC and ATEEC. The back-to-back region cos <p -1 is well described, while small 

discrepancies, at the level of 10%, are observed in the central region of the TEEC and for large cos(/> 
values. The description by Herwig-i-i- is poorer. 

The shape of the ATEEC is very similar to that observed at colliders, see Refs. [23-35], and well 
reproduced by Pythia 6 and Alpgen. 


7 Theoretical predictions and uncertainties 

In perturbative QCD (pQCD), according to the factorisation theorem [66], final-state observables can 
be expressed as a convolution of the partonic cross-sections, &, with the parton distribution functions. 
Thus, in this particular case, the TEEC distribution to leading order in the strong coupling constant, 
can be expressed as the three-jet, energy-weighted, differential cross-section in cos (p, normalised to the 
integrated two-jet cross-section. This can be schematically expressed as 

1 di: ^ 

o-d(cos0) ’ 


( 4 ) 













































Figure 3; The unfolded distributions for transverse energy-energy correlation (left) and its asymmetry (right) along 
with comparisons to MC expectations. The statistical uncertainties are shown with error bars, while the total 
experimental uncertainties are shown in a shaded band. 


where the transverse energy-energy weighted partonic cross-section, xt {i = 1,2) are the 

fractional longitudinal momenta of the initial-state partons, fa^{xi) and fa 2 {^ 2 ) are the PDF, and ® de¬ 
notes a convolution over the appropriate variables. The denominator of Eq. (4) is the integrated dijet 
cross-section used to normalise the TEEC. 


The pQCD NEO calculations of the TEEC and ATEEC distributions are performed using NEOJet-i-i- 
[2, 3] interfaced with the MSTW 2008 [67], CTIO [68], NNPDE 2.3 [69] and HERAPDE 1.5 [70] par- 
ton distribution functions at NNEO. Typically, 0(10^°) events are generated for these calculations. This 
involves the calculation of the 2 —> 3 partonic subprocesses at NEO accuracy and of the 2 —> 4 partonic 
subprocesses at tree level. In order to avoid the double collinear singularities appearing in the latter [38], 
the angular range is restricted to | cos 0| < 0.92. 


The renormalisation and factorisation scales, inherent in any pQCD calculation, are usually taken to 
reflect the typical transverse momentum of the process under investigation. Eor the TEEC and ATEEC 
calculations, they are taken to be 




Pti + PT2 
2 


(5) 


where pxi and pt 2 are the transverse momenta of the two leading jets. This is also the choice in Ref. [71]. 
The value of the strong coupling constant at a given scale is connected to as(mz) using the two-loop beta 
function [2, 3]. 


The NEO theoretical predictions are subsequently corrected for non-perturbative effects such as hadron- 
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eos^ 

TEEC 

Stat. 

JES 

JER 

Shower 

Pileup 

Unfolding 

(-1.00, -0.96) 

10.008 

0.008 

+0.033 

-0.034 

0.009 

0.037 

0.008 

0.008 

(-0.96, -0.92) 

0.8218 

0.0047 

+0.0044 

-0.0040 

0.0011 

0.0044 

0.0036 

0.0005 

(-0.92,-0.88) 

0.3848 

0.0029 

+0.0029 

-0.0026 

0.0006 

0.0028 

0.0028 

0.0002 

(-0.88,-0.84) 

0.2324 

0.0022 

+0.0024 

-0.0022 

0.0004 

0.0023 

0.0022 

0.0001 

(-0.84, -0.80) 

0.1612 

0.0017 

+0.0022 

-0.0022 

0.0003 

0.0022 

0.0018 

0.0002 

(-0.80, -0.72) 

0.1095 

0.0009 

+0.0020 

-0.0020 

0.0002 

0.0020 

0.0015 

0.0002 

(-0.72, -0.64) 

0.0767 

0.0008 

+0.0017 

-0.0017 

0.0001 

0.0017 

0.0012 

0.0001 

(-0.64, -0.56) 

0.0574 

0.0006 

+0.0015 

-0.0015 

0.0001 

0.0015 

0.0009 

0.0001 

(-0.56, -0.48) 

0.0472 

0.0005 

+0.0014 

-0.0014 

0.0001 

0.0014 

0.0005 

0.0001 

(-0.48, -0.36) 

0.0400 

0.0004 

+0.0012 

-0.0013 

0.0001 

0.0012 

0.0003 

0.0001 

(-0.36, -0.24) 

0.0329 

0.0004 

+0.0011 

-0.0012 

0.0001 

0.0010 

0.0001 

0.0001 

(-0.24,-0.12) 

0.0302 

0.0003 

+0.0010 

-0.0011 

0.0001 

0.0009 

0.0001 

0.0001 

(-0.12,0.00) 

0.0273 

0.0003 

+0.0009 

-0.0010 

- 

0.0008 

0.0001 

0.0001 

(0.00,0.12) 

0.0262 

0.0003 

+0.0009 

-0.0010 

- 

0.0008 

0.0001 

- 

(0.12,0.24) 

0.0264 

0.0003 

+0.0009 

-0.0010 

- 

0.0008 

0.0002 

- 

(0.24,0.36) 

0.0272 

0.0003 

+0.0009 

-0.0010 

0.0001 

0.0009 

0.0004 

- 

(0.36,0.48) 

0.0286 

0.0003 

+0.0010 

-0.0010 

0.0001 

0.0010 

0.0006 

- 

(0.48,0.56) 

0.0306 

0.0004 

+0.0011 

-0.0010 

0.0001 

0.0011 

0.0008 

- 

(0.56,0.64) 

0.0340 

0.0004 

+0.0012 

-0.0011 

0.0001 

0.0011 

0.0006 

- 

(0.64,0.72) 

0.0391 

0.0004 

+0.0014 

-0.0014 

0.0001 

0.0012 

0.0004 

0.0001 

(0.72,0.80) 

0.0487 

0.0004 

+0.0017 

-0.0018 

0.0001 

0.0013 

0.0002 

0.0001 

(0.80,0.84) 

0.0639 

0.0007 

+0.0024 

-0.0026 

0.0001 

0.0014 

0.0002 

0.0002 

(0.84,0.88) 

0.0780 

0.0008 

+0.0029 

-0.0032 

0.0002 

0.0014 

0.0002 

0.0004 

(0.88,0.92) 

0.0955 

0.0009 

+0.0031 

-0.0033 

0.0002 

0.0013 

0.0003 

0.0005 

(0.92,0.96) 

0.1025 

0.0009 

+0.0021 

-0.0022 

0.0001 

0.0009 

0.0003 

0.0004 

(0.96,1.00) 

11.448 

0.003 

+0.039 

-0.036 

0.006 

0.030 

0.008 

0.008 


Table 1: Values of the transverse energy-energy correlation function (TEEC). The statistical and systematic uncer¬ 
tainties, due to Jet Energy Scale and Resolution (JES and JER), shower modelling as well as pileup and unfolding, 
are shown in the subsequent columns. Uncertainties marked with a dash (-) are smaller than 0.00005. 


isation and the underlying event. This eorreetion is ealeulated using the leading-logarithm parton shower 
generators Pythia 6 and Herwig-i-i- interfaeed with different tunes. The full MC generator partiele-level 
predietions with these effeets switehed on are eompared with the parton-level predietions before had- 
ronisation and without UE effeets. From this eomparison a bin-by-bin eorreetion faetor is ealeulated as 
the ratio of the two predietions, whieh is then used to eorreet the NLOJet-i-i- output. They are found to 
deviate from unity by about 1% for both Pythia 6 and Herwig-i--i- for most of the | eos (p\ < 0.92 range. 

Three main theoretieal uneertainties are eonsidered for the analysis: those eorresponding to the renor- 
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cos (p 

ATEEC 

Stat. 

JES 

JER 

Shower 

Pileup 

Unfolding 

(-1.00, -0.96) 

-1.4406 

0.0083 

+0.0094 

-0.0066 

0.0144 

0.0078 

0.0010 

0.0001 

(-0.96, -0.92) 

0.7193 

0.0048 

+0.0002 

-0.0000 

0.0002 

0.0012 

0.0037 

0.0001 

(-0.92, -0.88) 

0.2893 

0.0030 

+0.0012 

-0.0008 

0.0012 

0.0022 

0.0028 

0.0003 

(-0.88,-0.84) 

0.1544 

0.0023 

+0.0009 

-0.0006 

0.0007 

0.0019 

0.0023 

0.0002 

(-0.84, -0.80) 

0.0973 

0.0019 

+0.0007 

-0.0005 

0.0004 

0.0015 

0.0020 

- 

(-0.80, -0.72) 

0.0608 

0.0010 

+0.0006 

-0.0006 

0.0002 

0.0010 

0.0016 

- 

(-0.72, -0.64) 

0.0376 

0.0009 

+0.0005 

-0.0006 

0.0001 

0.0007 

0.0011 

0.0001 

(-0.64, -0.56) 

0.0235 

0.0007 

+0.0004 

-0.0005 

0.0001 

0.0004 

0.0007 

0.0001 

(-0.56, -0.48) 

0.0165 

0.0007 

+0.0004 

-0.0005 

0.0001 

0.0003 

0.0005 

0.0001 

(-0.48, -0.36) 

0.0115 

0.0005 

+0.0003 

-0.0004 

0.0001 

0.0002 

0.0004 

0.0001 

(-0.36, -0.24) 

0.0057 

0.0004 

+0.0001 

-0.0002 

- 

0.0001 

0.0002 

0.0001 

(-0.24, -0.12) 

0.0038 

0.0004 

+0.0001 

-0.0001 

- 

0.0001 

0.0001 

0.0001 

(-0.12,0.00) 

0.0011 

0.0004 

- 

- 

- 

- 

- 


Table 2: Values of the asymmetry on the transverse energy-energy correlation function (ATEEC). The statistical 
and systematic uncertainties, due to Jet Energy Scale and Resolution (JES and JER), shower modelling as well as 
pileup and unfolding, are shown in the subsequent columns. Uncertainties marked with a dash (-) are smaller than 
0.00005. 


malisation and factorisation scale variations, those corresponding to the PDF, and those on the non- 
perturbative corrections. 

• Scale uncertainty: The ambiguity in the choice of the renormalisation and factorisation scales gives 
rise to a scale uncertainty. To estimate it, the scales jur and /ip are varied by a factor of two up and 
down, with the additional requirement that 0.5 < /tr/z/f ^ 2. From all those variations, the largest 
uncertainty is obtained when both /ir and jur are varied simultaneously by the same factor from the 
nominal scale. These two combinations are used to define the envelope of the scale uncertainty for 
both the TEEC and ATEEC. The size of the scale uncertainty is highly asymmetric and is at most 
about 8% for the TEEC distribution, and somewhat smaller for the ATEEC. 

• PDE uncertainty: The CTIO parton distribution functions provide 50 variations for the 25 fitted 
parameters at the 90% confidence level. Each of fhe 25 paramefers are varied up and down fol¬ 
lowing fhe CTIO recommendations in Ref. [68], and are combined for each bin of fhe TEEC and 
ATEEC disfribufions following fhe prescripfion given in Ref. [72]. The size of fhe PDE uncer- 
fainfy, once scaled af 68% confidence level, is abouf 1.5%. A similar procedure is used for fhe 
MSTW2008, NNPDE 2.3 and HERAPDE 1.5 parfon disfribufion functions. 

• Uncerfainfies in fhe non-perfurbafive corrections: The non-perfurbafive corrections (NPC) are cal- 
culafed using Pythia 6 inferfaced fo fhe AUET2B and AMBT2B funes [45, 46], as well as Her- 
wiG-i-t wifh fhe UE7000 fune [50]. Moreover, Pythia 8 inferfaced fo fhe 4C and AU2 funes is also 
used. An uncerfainfy is derived by considering, on a bin-by-bin basis, fhe maximum difference 
befween fhe nominal Pythia AUET2B and any ofher fune. Its size is below 1% for most of the 
angular range considered. 
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8 Determination of the strong coupling ors(mz) 


The evaluation of a^{mz) is made by minimising a;y^ function taking into account correlations between 
the systematic uncertainties using nuisance parameters Ak, one for each source of uncertainty. These 
nuisance parameters are normalised to zero mean and unit variance. The minimum of the function 
is found in a 66-dimensional space, one dimension corresponding to a^{mz) and the rest to the nuisance 
parameters associated with the experimental errors. The function to be minimised is defined as 


i 


(Xj - Fija^, 

Axf + Ar^ 



where the NLOJet-i-i- predictions are varied according to 

Fi(as, A) = iA,(o's) 

K k ) 


1 + 2 ] Ako-f 


( 6 ) 


(V) 


In these expressions, x, corresponds to the data points in each distribution (TEEC or ATEEC), and Ax; 
are their statistical uncertainties. At,- are the statistical errors on the NEOJet-i-i- predictions, while 
correspond to the k-th source of experimental uncertainty in the bin i. 

The functions i/^,(as) are analytical expressions parameterising the dependence of each observable (TEEC 
or ATEEC) on the strong coupling constant. They are obtained by fitting the predictions for each bin as 
a function of as(wJz)- This function is chosen to be a parabola, as the theoretical predictions account for 
terms quadratic in as- The quality of the fit to the NEO theoretical predictions is found to be excellent for 
each bin of the TEEC and ATEEC. The uncertainties from these fits are negligible. 

The theoretical uncertainties on the predictions are treated by varying the theoretical distributions by each 
independent source of uncertainty (scale, all independent PDE uncertainties and non-perturbative correc¬ 
tions) and repeating the fit using the modified theoretical input. 


The fit to the TEEC data exhibits shifts in a few nuisance parameters, which are always compatible with 
the ±lcr band. The results for the strong coupling constant obtained using different parameterisations 
of the PDE are summarised in Table 3, together with the experimental uncertainties and the values of 
X^INdof- 


PDE 

as(mz) value 

X^INdof 

MSTW 2008 

CTIO 

NNPDE 2.3 

HERAPDE 1.5 

0.1175 + 0.0010 (exp.) (scale) + 0.0006 (PDE) + 0.0002 (NPC) 

0.1173 + 0.0010 (exp.) (scale) + 0.0017 (PDE) + 0.0002 (NPC) 

0.1183 + 0.0010 (exp.) (scale) + 0.0009 (PDE) + 0.0002 (NPC) 

0.1167 + 0.0007 (exp.) (scale) (PDF) + 0.0001 (NPC) 

29.0/21 

22, A 111 

29.3/21 

22.1 111 


Table 3; Results for from fits to the TEEC function using different PDE sets, namely MSTW 2008, CTIO, 
NNPDE 2.3 and HERAPDE 1.5, together with experimental as well as theoretical uncertainties due to scale and 
PDE choices and non-perturbative corrections. 


The final value for the TEEC fits is chosen to be the one obtained using CTIO, since its PDE uncertainty 
is largest and serves as an envelope covering the variations with different PDE sets as shown in Table 3: 
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Figure 4: The unfolded distributions for transverse energy-energy correlation (left) and its asymmetry (right) com¬ 
pared with the results of a fit to pQCD NLO calculations including non-perturbative corrections. The green shaded 
band indicates the uncertainty on the theoretical predictions, which includes the sum in quadrature of uncertainties 
associated with scale, a^, PDF and NPC. The statistical uncertainties on the predictions are indicated by green error 
bars, appreciable only on the tail of the ATEEC. The solid error bars on the data points (in black) indicate the ex¬ 
perimental uncertainties taking into account the correlations between them. The fitted values of the strong coupling 
constant are a^\inz) - 0.1173 (TEEC) and af\mz) = 0.1195 (ATEEC). 


as(mz) - 0.1173 ± 0.0010 (exp.) (scale) ± 0.0017 (PDF) ± 0.0002 (NPC) . (8) 

The fit to the ATEEC data does not show any significant shift in the values of the nuisance parameters. In 
this case, the fit results in the values for the strong coupling constant which are summarised in Table 4. 


PDF 

ajjriz) value 

A^/Ndof 

MSTW 2008 

CTIO 

NNPDF 2.3 

HERAPDF 1.5 

0.1195 + 0.0017 (exp.) (scale) + 0.0006 (PDF) 

0.1195 + 0.0018 (exp.) (scale) + 0.0016 (PDF) 

0.1206 + 0.0018 (exp.) (scale) + 0.0009 (PDF) 

0.1182 + 0.0013 (exp.) ro.ooov 

12.7/10 

12.6/10 

12.2/10 

12.1/10 


Table 4; Results for Us from fits to the ATEEC function using different PDE sets, namely MSTW 2008, CTIO, 
NNPDE 2.3 and HERAPDE 1.5, together with experimental as well as theoretical uncertainties due to scale and 
PDE choices. The uncertainty due to the non-perturbative corrections is negligible. 

The final value for the ATEEC fit is also chosen to be the one obtained using the CTIO parton distribution 
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functions: 


a,{mz) = 0.1195 ± 0.0018 (exp.) (scale) ± 0.0016 (PDF). 


(9) 


The agreement between the fitted theoretical NLO predictions, including non-perturbative corrections, 
and the data is good as shown in Fig. 4 and indicated by the;y^ values given in Tables 3 and 4. Restricting 
the angular region in the fits to (-0.72,0.72), yield values of the strong coupling constant which vary 
within experimental uncertainties. The values of a^imz) found in this analysis are in agreement with the 
world average a^{mz) - 0.1185 ±0.0006 [73], as well as with other determinations of the strong coupling 
constant from the data collected at the LHC [71, 10, 74]. 

Calculations beyond NLO accuracy, which are already available for processes such as top-quark pair [75] 
or Higgs boson production [76], are needed for multi-jet production at LHC energies. They are expected 
to reduce the scale uncertainties, which are the limiting factor in this determination of the strong coupling 
constant. 


9 Summary 


First measurements of the TEEC and ATEEC functions are presented using 158 pb“^ of pp collision 
data at 7 TeV recorded by the ATEAS experiment at the EHC. Eor this purpose, multi-jet final states 
are selected requiring jets, reconstructed with the anti-k, algorithm and radius parameter R - 0.4, with 
Pt > 50 GeV and \r}\ < 2.5 and such that the scalar sum of the transverse momenta of the two leading jets 
is above 500 GeV. The TEEC and ATEEC data are fairly well described by Pythia 6 and Alpgen, while 
the Herwig±± MC simulation shows some discrepancies which can be as large as 30%. 

The TEEC and the ATEEC at the particle level are compared to perturbative QCD predictions at NEO 
accuracy. The renormalisation and factorisation scales are chosen to be (pxi + Pt 2)/2, ranging from 250 to 
1300 GeV and with an average value of 305 GeV. Through their construction, both the TEEC and ATEEC 
functions are less affected by experimental effects such as the jet energy scale and resolution or pileup than 
absolute cross-section measurements. Similarly, the PDE uncertainties in their theoretical predictions, as 
given by Eq. (4), cancel to a large extent. This renders these observables well suited to determine the 
strong coupling constant. The data for | cos (l)\ < 0.92 are fitted to the QCD predictions obtained with 
NEOJet±± to determine the value of the strong coupling constant. Eor the TEEC, which provides the 
experimentally more accurate determination, the result of the fit using the CTIO PDE yields 

a^(niz) = 0.1173 ± 0.0010 (exp.) (scale) ± 0.0017 (PDE) ± 0.0002 (NPC). (10) 

The present determination of Osimz) is limited by the uncertainties due to the choice of renormalisation 
and factorisation scales. 
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